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CN ■ Abstract 
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' Contaldi et al. [I] have suggested that an initial period of kinetic energy domination in 

. single field inflation may explain the lack of CMB power at large angular scales. We note 

that in this situation it is natural that there also be a spatial gradient in the initial value 
of the inflaton field, and that this can provide a spatial asymmetry in the observed CMB 
power spectrum, manifest at low values of t. We investigate the nature of this asymmetry 
C^i . and comment on its relation to possible anomalies at low £. 

6 

& 

in 

> 



1 E-mail: donoghue@physics .umass . edu 
E-mail: koushik@mppmu.mpg.de 
3 E-mail: andreas.ross@yale.edu 



1 Introduction 



The simplest models of inflation involve only a single scalar field, the inflaton [2]. With a 
suitably chosen potential, such a model can provide a simple explanation of the temperature 
fluctuations in the CMB at all angular scales [3]. There are however a number of possible 
anomalies in the CMB power spectrum that have attracted the attention of many researchers 
@], [!], [6], [7], [8], [9], [10], p], P2], [J3], [U]. (However see also [H].) None of these 
anomalies is by itself statistically compelling; however, taken together they provide a hint 
that these features may be significant. Much discussion of anomalies involves the power 
spectrum at low £, i.e. at large scales, where several anomalies indicate a possible spatial 
asymmetry in the power spectrum, most often roughly a north-south galactic coordinate 
asymmetry [4], 0, ©, [9], [UJ, [12], [JJ], [M]. 

The possibility that there may be an asymmetry in the observed CMB power spectrum 
was first raised by Eriksen et al [4] and Hansen et al [5] using the first year WMAP data. 
Their data analysis suggested a difference in power of roughly 20% for low i maximized in 
the direction of galactic coordinates (80°, 57°). Interestingly no effect was seen above i ~ 40. 
For example, the analysis of the power spectrum in the vicinity of the first acoustic peak 
[5j [16] showed no evidence of a spatial asymmetry. At low values of £, the cosmic variance 
provides an intrinsic scatter in the power spectrum data, so that even though the signal is 
rather large, the statistical significance of their result was below 3 standard deviations. 

In their three year data release [3], the WMAP team addressed the isotropy of the power 
spectrum, finding a small asymmetry in a direction consistent with Eriksen et al [4J. The 
method introduced by the WMAP team to investigate asymmetries in the CMB spectrum 
is to multiply an isotropic Gaussian CMB field by a large scale modulation function. They 
test both a dipole and a quadrupole modulation and find that the significance of the signal is 
not statistically compelling. Their analysis uses a pixel size of 7° which makes their analysis 
sensitive up to I ~ 20. The original Eriksen team has also revisited the WMAP 3 year data 
[9j using a statistical framework similar to the WMAP team's with a modulation function. 
They choose a higher resolution with a pixel size of 3.6° including multipoles up to I ~ 40 
in their analysis and confirm the asymmetry with a higher statistical significance than the 
WMAP team and in consistency with their previous analysis of the first year WMAP data. 

Hansen et al [UJ and Maino et al [12] explored two different approaches to extract the 
CMB spectrum where WMAP data itself is used for foreground removal, and both find an 
asymmetry of the power spectrum at largest scales consistent with previous analysis and 
with each other. The fact that the asymmetry does not vary when different foreground 
subtraction procedures are applied constitutes a strong argument against a galactic origin for 
the asymmetry. Moreover, the asymmetry was also found in COBE data [3] which indicates 
that systematics may not be the correct explanation for a large scale asymmetry in the CMB 
power spectrum. Rath et al [13] have also found the asymmetry in the WMAP 3 year data 
using statistical techniques different from the ones used in previous analyses. Recent analyses 
of WMAP 5 year data show a similar anisotropy of power between the two hemispheres, but 
with the asymmetry possibly reaching to higher multipoles [14j . Otherwise, the nature of 
the asymmetry and the maximum asymmetry direction remains almost the same, and it will 
be interesting to see what results from Planck will have to say about an asymmetry at small 
scales. 

These analyses provide motivations for the study of inflationary models that can generate 
a spatial asymmetry at low t while remaining isotropic at larger values of I. If these anomalies 
prove to be valid indicators of an asymmetry in the power spectrum, they can provide a direct 
probe of inflationary dynamics. Significant work that attempts to find a solution to these 
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anomalies has already appeared [T7], [18], [19], [20] in the literature. 

In this paper we discuss a simple situation that could lead to a spatial asymmetry in the 
CMB power spectrum at low values of I within single field inflation. This involves an initial 
period of fast-roll expansion driven by the inflaton kinetic energy. The possibility of such an 
initial fast-roll period has been proposed by Contaldi, Peloso, Kofman and Linde (CPKL) 
[I] as a mechanism to explain the lack of CMB power at low I. This mechanism provides a 
suppression of the spectrum of primordial perturbations and thus of the CMB at large scales, 
and it has also been worked on by others [21], 

We will argue that in situations where the initial kinetic energy is significant in comparison 
to the potential energy, we should also expect the presence of a spatial gradient in the initial 
conditions of the inflaton field. We will show that even a surprisingly small value of an initial 
gradient - of order a few percent - will leave an observable spatial asymmetry in the CMB 
power spectrum at low L Essentially, the power suppression in the fast-roll model occurs at 
scales that depend sensitively on the initial magnitude of the scalar field in the frame where 
the kinetic energy is uniform and isotropic. This leads to a characteristic pattern for the 
spatial dependence of the power spectrum. While we will provide a brief discussion of two- 
field models below, we here focus on the single field fast-roll option because of its simplicity 
and predictive power. 



2 Kinetic energy and spatial asymmetries 

Inflation provides an explanation for the isotropy and homogeneity of the present universe. 
Rather than having to postulate extremely smooth initial conditions for the early universe, 
a long period of inflation will take non-smooth initial conditions and still lead to a highly 
isotropic and homogeneous observable universe today. However, if the number of e-foldings of 
inflation is just barely the minimal number, about 60, the initial conditions could be relevant 
and could modify the first few e-foldings. 

The CPKL mechanism [I] postulates an initial period of kinetic energy dominance which 
then rapidly evolves into the standard slow-roll paradigm where the potential energy dom- 
inates and the universe inflates. If the slow-roll phase is many e-foldings longer than the 
minimum number of e-folds, the effects of the initial kinetic phase will be unobservable since 
the scales associated to its effects will be stretched far beyond our observable horizon. How- 
ever, if the slow-roll phase is close to the minimum number of e-folds, then that initial kinetic 
phase will modify the first few e-foldings that generate the CMB power spectrum on largest 
scales, i.e. for small values of I. 

In a universe dominated by a uniform scalar field, the equations of motion are the Klein- 
Gordon equation on a FRW background 

j> + 3H<j> + ^-{<j)) = (1) 



and the Friedman equation^ 

H 2 ={-\ -^-tIt-MmI. (2) 



d\ 2 8vr 8vr 1 12 



aj 3 r 3 \ 2 
and the equation of state parameter is 



p \ - v(<j>) 

w = - = -f 2 • (3) 



x We use units where M P = G~ 1/2 = 1. 
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The CPKL assumption is an initial condition in which the inflaton velocity, eft, is non-zero 

and the kinetic energy term is dominant over the potential energy, ^- 3> V{4>). During this 
initial phase w ~ 1 and the expansion of the universe will be decelerating similar to a matter 
dominated universe rather than a deSitter expansion. In this phase the kinetic energy rapidly 
decreases, until eventually the potential energy dominates and we enter the usual slow-roll 
phase. There is a short transitional phase when the potential energy already dominates 
and the universe inflates, but the inflaton velocity has not yet settled to its slow-roll value 

<PSR ~ 2y^ d<t>l V ■ 

The initial conditions involve specifying both <j> and <p which have to be chosen appropri- 
ately to obtain only about 60 e-folds of inflation so that the effects of the initial fast-roll stage 
are observable. CPKL then show that the quantum fluctuations are suppressed during the 
onset of inflation when the inflaton is fast-rolling - this will be reviewed in the next section. 
The picture that emerges then involves suppressed fluctuations at early times followed by 
standard slow-roll behavior. Since the earliest times correspond to the largest scales, the low 
I multipoles are suppressed while the higher ones are standard. 

Our extension of CPKL comes from the observation that the initial conditions in <f> and <f> 
need not be the same at all positions in space. If they are close to uniform, one can expand 
the values in a multipole expansion. The first deviation from uniformity would consist of a 
gradient in the initial conditions across the initial patch. We will consider only such leading 
linear deviations from uniformity in this paper. 

CPKL invoke a uniform initial condition in (j). Actually this is not a separate assumption. 
Because the value of <j) is changing with time, one can always choose a time-slicing such that 
4> is uniform across the initial time slice. That is, if there is a gradient in the initial condition 
for cf) using one definition of the initial time slice, one can change to another definition such 
that this variable is uniform. However, in this frame there is no a priori reason for the initial 
value of the magnitude of (j) itself to be uniform. A mechanism that can produce a temporal 
variation in (ft can in principle also produce a spatial gradient in the field. We could equally 
well define a different time frame in which the initial condition of <f> is uniform, but in this 
frame we would in general not expect that the initial value of (f> is constant in space. It is an 
extra assumption to assume that the initial conditions for and <j> are spatially uniform in 
the same frame. 

We can obtain an invariant description of the evolution of the inflaton by considering 
trajectories in the <f> — (f> plane, displayed in Fig. Q] for a chaotic inflation potential V((j)) = 
^m 2 (j) 2 . There are three phases visible in this plot: Initially, kinetic energy dominates and 
due to the rapid decrease in the kinetic energy, the trajectory runs quickly to the slow-roll 
attractor line where <j> ~ (j)sR = const. All initial trajectories are attracted to this line - this 
is one of the key features of inflation. Finally at the end of inflation, there is inflaton decay 
and reheating. However, different starting points lead to differing amounts of the slow-roll 
phase - these are shown as different initial trajectories. The number of e-folds of slow-roll 
inflation is 

N * 2ir</> 2 SR>i (4) 

where (j)sR,i is the field value when the phase space trajectory hits the slow-roll attractor. 

If we start in a frame with a uniform value of 4>, which would be represented by a vertical 
slice through the initial trajectories, one needs significantly different initial values of <fi in 
order to end up at different points on the slow-roll attractor. However, if we use a frame 
with uniform initial values of (j), which would be a horizontal slice through the trajectories, 
only a small difference in the initial values of <j> are required to produce distinct trajectories 
with differing amounts of slow-roll behavior and thus differing amounts of inflation. We will 
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Figure 1: Phase space diagram of inflationary background solutions for two different values of the 
initial scalar field with both initially in the kinetic energy domination regime. The dashed curve shows 
the slow-roll attractor. 

see that gradients of order only a few percent are needed for observable effects. Because the 
potential energy is subdominant at the initial time, this is only a very small gradient in the 
initial energy density. Therefore, an initial slice with constant (j> is the better choice since we 
want to work with a FRW metric that requires a homogeneous and isotropic energy density. 

If the inflaton field is not uniform, the Klein-Gordon equation, Eq. (pQ), contains an 
additional term proportional to V 2 (p. For a gradient in the field, 4>(x) = A + Bx where A 
and B are constants, V 2 </> = so that we can neglect this term and still work with Eq. (JTJ) . 
The gradient however contributes a term to the energy density of the form (V0) 2 /2, which 
for a linear gradient in the field (j) amounts to a constant throughout space. There is also 
an anisotropic contribution to the pressure proportional to (Vj</>) 2 . However, the gradient's 
contribution to the energy density and pressure is always subdominant for the small amounts 
of gradients we require so that we can neglect it. Thus, even in the presence of a small 
spatial gradient, the inflaton field evolves independently at each spatial position. That is, in 
our approximation the trajectories displayed in Fig. 1 are not modified by the presence of a 
spatial gradient. 

In our extension with a gradient in the initial conditions, two points on opposite sides 
of the universe which started with different initial values of <p will have undergone different 
amounts of inflation so that the large scale suppression features in the power spectrum asso- 
ciated to an initial fast-roll stage will appear at different scales today. A gradient in e-folds 
of inflation is the leading effect in our model stretching both the cutoff scale and space by a 
different amount in different parts of the universe. Certainly, it is not strictly correct to use 
a FRW background metric and the standard formulas for a uniform and isotropic cosmology. 
However, the expansion is uniform and isotropic both during the fast roll phase when the 
kinetic energy is dominant, and later during the slow roll phase. The inhomogeneity only 
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effects the universe during the very short transition region between fast-roll and slow-roll of 
the inflaton. We feel that our treatment captures the leading effects of an initial gradient 
without the need to solve exactly the evolution through the short transition region. 



3 CMB fluctuations at low multipoles 

An initial regime of kinetic energy dominance of the inflaton before reaching the slow-roll 
attractor modifies the spectrum of quantum fluctuations as noted by CPKL. For wavelengths 
much smaller than the Hubble scale at the onset of the accelerated expansion, H^ 1 , the 
spectrum is the same as the slow-roll spectrum because the small wavelength perturbations 
are insensitive to the overall background expansion of the Universe at this time. However, 
for wavelengths comparable to H^ 1 the spectrum is altered and it exhibits a suppression 
of power for larger wavelengths. Intuitively, this suppression can be understood from the 
relation Sp/p ~ H 2 /<fi: if the inflaton rolls faster initially, the spectrum will be suppressed at 
larger scales. 

Following CPKL, we use a chaotic inflation model with potential V(4>) = ^m 2 (j) 2 with 
m = 4 x 10~ 7 . As initial conditions, we choose 4>i = 3.59 and 4>i = —35.9m. This particular 
choice of initial conditions gives us about 60 e- folds of inflation. The initial kinetic energy 
is 100 times larger than the initial potential energy so that we start out well in the kinetic 
energy domination regime. 

We use the gauge invariant formalism for cosmological perturbations [221 [23] and the 
equations of motion for the perturbation variables read 



5<j>l + 2H6(f> k +( , ^\ 8<j> k + V M a 2 8<l) k = 40^ - 2V 4> a 2 $> k (5) 

4:ir<f>' 5<j> k = $' fc + ft$fe (6) 

where primes denote differentiation with respect to conformal time r/ and TL = a /a. We solve 
these perturbation equations using numerical mode by mode integration. As the system is 
initially in the kinetic energy domination regime, the vacuum is chosen to be the approximate 
solution in the kinetic energy domination stage at an initial time rji pQ 



v k (vi) = sfi + m^i (km + ^ ) , in 

where the Mukhanov- Sasaki variable is 

v k = a I 5cj) k + I . (8) 




We verified that a different choice of vacuum, the Hamiltonian diagonalization vacuum given 
by [22] 

/ //\-l/4 / „ \ 1/4 

v k (m) = [k 2 - Z —\ , v' k ( Vi ) = i Ik 2 - Z —\ withz^^- (9) 

does not change the resulting power spectrum significantly. 

In Fig. [2] we show the power spectrum of the gravitational potential after the end of 
inflation. The cutoff of the spectrum around k ~ cloHq is the most important feature of 
the spectrum, but some oscillations at the transition to a pure slow-roll spectrum are visible 
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Figure 2: Power spectrum of the gravitational potential $ at the end of inflation in the CPKL model 
with an initial fast-roll stage. 

as well. The position of the cutoff in the spectrum, i.e. the scale associated to the cutoff, 
actually depends on the complete expansion history of the universe since the onset of inflation 
until today, which is not exactly known. 

Instead of computing the spectrum numerically CPKL have also shown that an approxi- 
mate analytic solution with an instantaneous transition between the regimes of kinetic energy 
domination and slow-roll inflation reproduces the exact spectrum extremely well. This shows 
that the spectrum does not depend much on the details of the innaton potential. We have 
checked this by generating the corresponding spectra using 4 and <f> e potentials, and we 
found shapes of the resulting power spectra similar to the one shown in Fig. [2] with only 
small changes in the width of the cutoff region of the order of 10%. 

Since the main feature of the spectrum is the cutoff that yields suppression for largest 
scales, CPKL have also introduced a simple parametrization of a spectrum with an exponen- 
tial cutoff 

P*(k) ~k 3 \$ k \ 2 =A s (l-exv[-(k/k c ) a )}k n °- 1 (10) 

with a = 3.35 useful simplified model. 

The spectrum of CMB fluctuations can be computed from the primordial power spectrum 
P$(fc). A variety of astrophysical processes make this a complicated and highly numerical 
task which is usually performed with numerical codes \24 \ 125 1 [26] . Since the power spectrum 
in wavenumber is not modified at large k, the CMB spectrum will be the same as the standard 
slow-roll predictions at high multipoles £. At low i the leading effects are the Sachs Wolfe 
effect and the late time integrated Sachs Wolfe effect. For the calculation of the CMB 
spectrum we have modified the numerical code CAMB [26J to include the CPKL primordial 
power spectrum as the initial spectrum. The resulting CMB spectrum is shown in Fig. [3l 
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Figure 3: TT power spectrum of the CMB in the CPKL model. 

One of the main results is the fact that the features in the primordial spectrum get somewhat 
smoothed in the CMB spectrum but the shape is generally not altered very significantly. We 
see the power suppression at low multipoles due to the initial fast-roll stage of the inflaton 
field, and at high multipoles the spectrum matches the flat slow-roll spectrum. The effect of 
the wiggles around the cutoff is also evident in the CMB spectrum. However, we note that 
unlike the primordial spectrum of Fig. [2j the CMB spectrum for multipoles smaller than 
the scale of the cutoff does not approach zero as quickly, but seems to have an offset which 
depends on the scale of the cutoff, see also Fig. HI This offset arises due to the late time 
integrated Sachs- Wolfe effect which gives contributions to Ce from a wider range of scales 
than the Sachs- Wolfe effect does. In the presence of a cutoff in the spectrum at large scales, 
the late time integrated Sachs- Wolfe effect is therefore more important than the Sachs- Wolfe 
effect at very low £ [27] . 

Since the details of inflation and reheating are not pinned down precisely, the expansion 
history of the universe is not known exactly and the number of e-folds N during inflation 
enters as an adjustable parameter, where a minimum of about 60 is needed to solve the 
flatness and horizon problems. Changing N in our case will shift the position of the feature 
in the spectrum. For N 3> 60, the feature is at scales much larger than our present horizon 
and the observable CMB spectrum is indistinguishable from a slow-roll spectrum. If two 
different parts of the universe underwent different amounts of inflation N after the initial 
fast-roll stage due to different initial conditions fa and fa, their spectra will match at small 
scales or large £ as seen in Fig. HJ The power spectrum in the part of the universe that 
inflated more e-folds will have the scale associated to the feature in the spectrum stretched 
more so that the feature appears at smaller multipoles £. There is in addition a geometric 
effect due to a temporary asymmetric expansion, which we explore in the next section. The 
CPKL mechanism is reasonably predictive in that there is only one free parameter, the scale 
of the feature in the spectrum, which is determined by the number of e-folds N of inflation 
and thus by the relative values of the initial conditions for <p an d <fi - see Fig. [Hand Eq. (jlj). 
Observationally, this is manifest as the value of £ below which there is a suppression in the 
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Figure 4: CMB spectra for different numbers of e-folds of inflation. The solid line comes from N 3> 60 
so that the cutoff feature due to an initial fast-roll stage is stretched beyond our observable universe 
and thus it corresponds to a pure slow-roll spectrum. 

CMB power spectrum. 



4 Gradient phenomenology and hemisphere analysis 

With a spatial gradient in the initial value of 4> in the frame where <ft is uniform, parts 
with different initial values undergo different amounts of inflation. Once all portions of the 
observable universe are in the slow-roll phase, the standard inflationary description holds 
with CMB fluctuations at these scales and subsequent structure formation being spatially 
uniform. The effect of the initial gradient only modifies the transition region between kinetic 
domination and slow-roll. This transition yields the cutoff in the spectrum at the scale k 
that will be stretched by different factors in different parts of the sky due to the presence of 
the gradient. 

Let us consider two different parts of the universe with slightly different initial field values 
<f>i t i and <j>i 2 = <f>n + A(f>. This also gives different amounts of inflation in the two parts, Ni 
and N 2 = Ni + AN where the difference in e-folds scales proportionally to A<p 



AN ~ A(f>. 



Normalizing the scale factor at the initial time, a different amount of e-folds of inflation gives 



a 2 /ai 



„A7V 



after inflation so that the relationship between coordinates and wavenumbers 



in the early universe and physical scales today then depends exponentially on the difference 
in e-folds. 

Our first task is to understand how a linear gradient in the scalar field translates to 
variations in the physical scales on the surface of last scatter. Let the coordinate C describe 
the direction along which the initial value of the inflaton field has a linear gradient in the 
frame where (f> is constant, i.e. 

MC) = A + B(. (12) 
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A small patch of the initial volume with thickness dC, will innate to a patch of the sky 
today with a thickness dz in physical coordinates, where z measures distance along the 
same direction today. Because different patches will have different amounts of inflation, with 
0.2 /di = e^ N and AN ~ A(f>, the thicknesses will be related by 

dz = pe K dC- (13) 

Here p is the rescaling of coordinates due to expansion that would have happened without 
the gradient in the scalar field (with fa = fa{(, = 0) = A): p ~ a(to)/a(ti). The parameter b 
depends on the magnitude of the gradient ^ in the inflaton field, with b = if there is no 
gradient in the field. Eq. (|13h shows the geometric effect of the gradient on the expansion 
of space: An initial patch located at a higher fa than a second patch and separated by a 
distance A£ will inflate AN = bA( e-folds more than the second one, so that the first patch 
will today be larger in the z-direction by an exponential factor of e^ N . 

The relation in Eq. (|13|) can be integrated to relate the initial coordinates along the 
gradient to coordinates today, 

z = £(e* - 1) (14) 

or 

bz~ 



bC = In 



1 + — 

P 



\n[l + b'z\. (15) 



Correspondingly we can relate a feature in the initial wavenumber spectrum, for example the 
location of the start of the cutoff in wavenumber, to the physical scales today. A cutoff that 
appears at k c in the original spectrum would appear at scales 

k = \ ee foe"* (16) 

today. Rewriting this in terms of the present physical coordinates implies that the feature 
appears at 

1 + a z 

It is apparent that due to the geometric effect of different amounts of expansion of different 
parts of the universe an initial gradient does not yield a simple gradient today! Finally, 
since we are interested in the CMB radiation coming from the surface of last scatter, we are 
interested in how this feature varies with direction. If the surface of last scatter is a distance 
i?* away, and we use z = cos 6 and define a = b'R*, we have the break in the spectrum 
occurring at 

k = -T J ^- Q - (18) 
1 + a cos 9 K ' 

The result is that an initial gradient in the scalar field, or equivalently in the start of the 
slow-roll phase, leads to the break in the CMB power spectrum occurring at a wavenumber 
that depends on the position in the sky by the above relation. 

It is tempting to translate Eq. (]18p into a relation for the cutoff position in terms of 
angles or multipoles in the CMB spectrum and use it to find the CMB modulation function 
for this model with a gradient. This is quite easy to do using the approximate primordial 
CPKL spectrum from Eq. (|10p and only taking into account the Sachs- Wolfe effect which 
results in a Ce spectrum of the same form as Eq. (flOl) (see the second toy model analysis in 
the Appendix). However as we have pointed out in the last section, the late time integrated 
Sachs- Wolfe effect is important at low £, and its inclusion unfortunately does not seem to 
result in a simple analytic modulation function. 
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Figure 5: Sketch of the measured CMB power spectra (mimicking the right panel of Fig. 7 in 
Ref. [12]) for the northern hemisphere (connected by dashed lines) and for the southern hemisphere 
(connected by solid lines), and the WMAP-3 best fit ACDM spectrum C° (smooth blue curve) [2"5] . 

A more rigorous test of the model would be to calculate the general correlations (a}, m ,a£ m ) 
including off-diagonal terms which vanish in the homogeneous and isotropic limit and compare 
these correlations to CMB data, as proposed for a different anisotropic model in [T7J. Such an 
analysis would go beyond the scope of this paper. Instead we will perform a simple first test 
of our model making use of the existing data in which the CMB power spectrum is extracted 
separately from two hemispheres of the sky. 

Data in hemisphere form has been reported by \TT\ [T2] and we sketch the results of 
the analysis of the WMAP 3-year data of [12] in Fig. [5] The two measured CMB spectra 
displayed are extracted from the two hemispheres with respect to the maximum asymmetry 
axis pointing to the north pole (9,4>) = (80°, 57°) in galactic coordinates, and the power 
spectrum for t < 40 obtained from the northern hemisphere exhibits a lack of power compared 
to the power spectrum of the southern hemisphere. 

In order to perform a first test of our model with an initial fast-roll stage and a spatial 
gradient in the initial field value, we approximate our results by performing a hemisphere 
averaging of the CMB spectrum and compare our predictions to the measured data in Fig. 
[5] For that we orient our initial gradient in the direction of the maximum asymmetry axis 
observed. Furthermore, we identify the point of the gradient with the lowest field value 
(where we expect the least amount of inflation and thus a cutoff feature in the spectrum 
present at smaller scales or higher I than the rest of the gradient) as the north pole. That 
then has the potential to yield a spectrum close to the observed one with a suppression in 
the northern hemisphere. 

Due to the fact that parts of the universe with different initial field values undergo different 
amounts of inflation, we first have to clarify what initial conditions lead to two equally sized 
hemispheres today. We start with the northern hemisphere (N) which contains the point with 
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Figure 6: The gradient in <pi shown in the original coordinates £ which are transformed into coordi- 
nates z in the present universe in which the surface of last scattering is shown. 



the lowest initial field value of the gradient, i.e. the lower part of the gradient, and choose 
a certain amount of initial gradient to give the northern hemisphere today. This setup is 
sketched in Fig. [6] in the coordinate Q and in today's physical coordinate z which are aligned 
with the gradient and the maximum asymmetry axis, and we normalize the coordinates such 
that the point £ = becomes z = today. If at Q = we expect N = Nq e-folds, at the north 
pole, i.e. the lowest point of the gradient Q = Clow, the amount of inflation is N = Nq — ANi ow . 
Since the number of e-folds scales linearly with the initial field value, the number of e-folds 
as a function of the spatial coordinate C is 

N(() = N -AN low -^. (19) 

Now we are interested in the size of this hemisphere today which follows from integrating 
dz = e^e N ^d^ where AT is the number of e-folds of the expansion of the universe since the 
end of inflation: 

zi ow = e M e N % low —— (20) 

which gives us the radius of the first hemisphere, i.e. the part of the gradient that inflated 
the least. As we want to construct two hemispheres of equal size, we require the size of the 
southern hemisphere (S) which contains the highest part of the gradient to be equal to the 
size of the northern hemisphere, 

Zhigh — %low \~-) 

With the relation between the coordinates £ an d z, we can in turn determine 

log 2 — e ~ ANlow 

Chigh = ~Clow a~\t (22) 

AJ-^low 



and find the number of e-folds for the point with most inflation from Eq. (|19p to be 

A^(C^)=iVo + log(2- e - AiV ^). (23) 

This result is interesting because it means that no matter how large a gradient there is in the 
northern hemisphere, the gradient in the southern hemisphere will at most give a difference 
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of log 2 ~ 0.69 e-folds of inflation within the southern hemisphere. This illustrates that parts 
of the universe with lower initial field values expand exponentially less than parts with higher 
field values so that the parts with highest initial field values dominate the universe today and 
parts with lowest initial field values comprise only a tiny fraction of the universe. 

For a power spectrum with a feature that varies as a function of position £ or z, we have to 
average the power spectrum over the two hemispheres in order to compare with the measured 
spectrum sketched in Figj5j Our approximation to this averaging is to imagine to cut each 
hemisphere into slices. If we use slices of equal thickness Az in today's coordinates, we take 
the average over all slices of the hemisphere by weighing each slice by its number of measured 
points on the surface of last scattering. The number of measured points is assumed to be 
proportional to the surface area of a slice, and since the surface area of a spherical segment 
only depends on its thickness Az and the radius of the sphere, each slice of equal thickness 
has the same weight in the average over the hemisphere. 

In the Appendix we present a simple toy model with a step function cutoff as the pri- 
mordial fluctuation spectrum. It is illustrative since most parts of the calculation can be 
performed analytically. Our analysis here however is performed completely numerically. We 
approximate our gradient by a series of equidistant steps and calculate the Ci spectra for 
each of the steps with a constant fa, i.e. we use equidistant (in coordinate £) slices of the 
gradient in which we approximate the initial field as constant. When averaging over a hemi- 
sphere, we have to average over all steps contained in the hemisphere. The geometric effect 
of different amounts of expansion of the different slices is accounted for by the weight factor 
exp ( —ANi ow j^— ) for each step, where we take the position C of the slice as its center value. 

\ slow J 

The effect of our mechanism on the averaged CMB spectra of the two hemispheres is best 
illustrated by defining relative suppression functions 

if = Cf/Cf* (24) 

with D = N,S as the ratio of the CMB spectra for the two hemispheres and the CAMB 
spectrum we obtain for pure slow-roll inflation. 

The two free parameters then in this hemisphere averaging procedure are the initial field 
values 4>i(C = 0) and fa(( = Clow), which then determine fa((, = Chigh) as well as Nq, ANi ow 
and AN^igh- In order to find a good fit to the hemisphere data in Fig. [51 we first generate 
a large array of slice spectra using CAMB. Here, we use slices which have a thickness A£ 
corresponding to AN s i ice = 0.1, i.e. the difference in number of e-folds between two adjacent 
slices is 0.1. In order to achieve the best fit to the data points for the two hemispheres, we 
vary the number of slices for the hemispheres (mostly the number of slices for the northern 
hemisphere since ANhigh ~ 0.7 for almost all of parameter space) and we vary the location of 
the border between the two hemispheres within our array of slices. That then corresponds to 
varying fa(( = 0) and fa(( = (i ow )- For each possibility, we obtain the suppression functions 
Fg and Ff for both hemispheres. These are then multiplied by the best fit ACDM curve 
C® for the WMAP 3-year data (see the smooth curve in Fig. [5]) to yield our mechanism's 
prediction for each hemisphere 

c Theory,D =F D C (25) 

with D = N, S. Here, (jJ heor y^ D — qD jf a rj parameters in CAMB are selected exactly as 
for the WMAP 3-year best fit. Since we use the slow-roll inflation spectrum in CAMB to 
obtain Cf R instead of a parameterized spectrum with a different spectral index, there could 
be small deviations due to the normalization of the CMB spectra, and we correct for these 
by using Eq. (|25p . The optimal configuration is then found from a \ 2 fit to the data points 
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Figure 7: Suppression functions resulting from averaged CPKL spectra for the northern (dashed) 
and southern (solid) hemisphere for ANi ow = 3.6. 



in Fig. [5] calculating 

x 2 = E E — "/ w) 

D=N,S I \°l ) 

where the variance is taken as the cosmic variance of the theoretical spectrum of the hemi- 
sphere, of = ^2/(2^ + 1) cJ heory ' D . 

The best-fit configuration has ANi ow = 3.6 and AN^^ = 0.7 which corresponds to less 
than 4% of variation in fa over both hemispheres together. In Fig. [7] we display the resulting 
suppression functions FP for both hemispheres. The spectrum of the southern hemisphere for 
the best fit is close to a slow-roll spectrum with little suppression only at very largest scales. 
For the northern hemisphere, the spectrum shows a suppression which is most significant at 
I < 10. Due to the large range of e-folds AN[ ow = 3.6, the averaged spectrum of the northern 
hemisphere does not exhibit any peak because the peaks of the individual spectra of the slices 
average out. 

When compared to the uniform and isotropic ACDM slow-roll CMB spectrum, we find 
that x 2 decreases, i.e. improves, by 3.1 for our model in the optimal configuration. Moreover, 
the prediction for the quadrupole £ = 2 for the full sky sphere is reduced by 45% in our 
model compared to the ACDM best fit prediction. However, this does not indicate that a 
significant preference for our fast-roll with gradient model arises: even though the decrease 
of Ax 2 = -3.1 improves the x 2 ■, t ne model introduces four new parameters when compared 
to the ACDM model - two angles for the direction of the gradient, the amount of gradient 
used AN[ ow and a scale fco at which the cutoff appears in a uniform spectrum of one of the 
slices. 

As one can see from Fig. [71 a significant suppression only arises for the lowest mutipoles 
I < 10 — 20. The reason for that is the geometric effect of the gradient which causes the 
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regions of the gradient with higher initial field values to dominate the universe, so that the 
resulting averaged hemisphere spectra still exhibit relatively steep cutoffs which are relatively 
close to each other. For the suppression functions of the hemispheres for the best fit shown in 
Fig. [7] the cutoff in the northern hemisphere's spectrum is roughly at £ ~ 2 whereas for the 
southern hemisphere it is at £ ~ 7. That makes it hard to improve the x 2 more significantly, 
since the hemisphere data sketched in Fig. [5] can be roughly characterized as a constant 
suppression of order 20% up to £ ~ 40. Moreover, the errors dominated by cosmic variance 
are reduced for suppressed theoretical spectra so that scattered data points can easily spoil 
the x 2 - 

In the Appendix we additionally discuss two toy models. The first one uses a step- 
function primordial spectrum, which allows us to perform parts of the hemisphere averaging 
analytically without having to approximate the gradient as steps. As a second toy model we 
consider a primordial spectrum of the same form as Eq. (jlOh but with a significantly softer 
cutoff than the CKLP spectrum. Such a more gradual cutoff could potentially achieve an 
asymmetry reaching to higher values of £. 

5 Conclusions 

The data appears to suggest a spatial modulation in the CMB power spectrum. However, this 
asymmetry is unusual in that it may only be present at largest scales or multipoles £ < 40. If 
inflation lasts only about 60 e- folds, this can occur if at the earliest stages of inflation there is 
an asymmetry which however disappears at later times of inflation. Provided the later stages 
of inflation are governed by a single slowly rolling field, the power spectrum at high £ and 
the universe today will be homogeneous and isotropic. However the quantum fluctuations at 
the earliest times, and hence the power spectrum at low multipoles, will show evidence of the 
initial lack of isotropy. 

We have provided a model of how this situation could have developed within the frame- 
work of single field inflation. It is based on the CPKL mechanism for large scale power 
suppression due to an initial fast-roll stage paired with a gradient in the initial values of the 
inflaton field. A gradient in the initial conditions is manifest as differing numbers of e-foldings 
of slow-roll behavior on different sides of the universe. This amounts to shifts in the £ values 
at which the suppression of the CMB power at large angular scales occurs. This model has a 
distinctive pattern for the generation of spatial asymmetries which makes it very predictive, 
and apart from the asymmetry axis, there are only two additional parameters characterizing 
the position of the cutoff and the amount of gradient within the observable universe. 

Our simple analysis where we approximated the initial gradient by a series of steps and 
averaged these over two hemispheres provided a first comparison to WMAP data extracted 
from different hemispheres. While our model reproduces the qualitative features of the data 
and improves the fit by A% 2 = —3.1, we cannot claim a preference for our model since it 
involves the introduction of in total four new parameters. Because of the geometric effect 
on the expansion of space resulting from differing amounts of inflation due to the initial 
gradient, a much better fit in a hemisphere analysis is hard to achieve with our simple single 
field model with an initial fast-roll stage and a gradient. Nevertheless, one could improve on 
our analysis by performing a full covariance analysis of the model including all anisotropic 
degrees of freedom, such as studied in [17j . In this way, one would use the information 
contained in the general correlations (fl^ m 'fl£m) ; including off-diagonal terms which vanish 
in the homogeneous and isotropic limit. 

It is also possible to generate spatial asymmetries in the power spectrum through the use 
of two or more fields. In this case, it is possible to assume that one field has a spatial gradient 
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in the frame in which the other field is uniform. Here there is a lot more freedom. With two 
fields, each field has its own potential with the possibility of a cross coupling between the 
fields. These possibilities enable one to modify the shape of the inflaton potential in different 
ways either enhancing or suppressing the rolling of the inflaton field and also modifying 
the amount of quantum fluctuations. The resulting possibilities for generating asymmetries 
deserve further study. Such models are more flexible, but are inherently less predictive than 
the single field model studied in this paper. However, at least some of these theories could 
have a similar phenomenology. By providing an initial faster evolution, one can suppress 
the curvature fluctuations providing a cutoff in the fluctuation spectrum, and a gradient in 
this initial condition would survive for only a few e-foldings. With the additional freedom in 
multi field models, one may be able to find a model where the geometric effect which limits 
our single field model's capability to obtain a significant asymmetry between hemispheres 
extending beyond I ~ 10 — 20 is absent or less restrictive. 

We note also that our model, and possible generalizations, have the potential to impact 
the analyses of other anomalies. For example, the cold spot uncovered in [8] occurs along 
the same axis as the power spectrum modulation. In our model the primary effect is the lack 
of power in some direction. Our hemisphere fit suggested that one hemisphere is close to 
the pure slow-roll while the other hemisphere shows the suppression. This could be similar 
in effect to the existence of a cold spot. Likewise, it is possible that this mechanism can 
modify the analysis of the unusual quadrupole-octopole alignments. When attempting a 
partial wave decomposition, an overall dipole shifts the underlying power from one multipole 
to neighboring ones with A£ = ±1. Our modulation in the spectrum is not exactly a dipole, 
but could nevertheless shift the apparent power from one multipole to nearby ones. 

It would be interesting to see if the addition of this form of spatial gradient to the studies 
that fit the CMB spectrum can provide further understanding of the proposed anomalies that 
are discussed in the literature. For that a more involved comparison to data than ours is 
needed. If the mechanism was successful it would add to our understanding of inflation and 
would increase the confidence in the inflationary paradigm. 
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Appendix 

In the Appendix, we study two toy models. The one considered first with a step function cutoff 
spectrum is instructive since the hemisphere averaging for the primordial power spectrum can 
be obtained analytically. For the second toy model we use an exponential cutoff spectrum 
which falls off more gradually than the CPKL spectrum. 

We construct the theta function toy model keeping the main aspects of the gradient fast- 
roll model, a large scale cutoff in the power spectrum whose scale varies within the universe 
due to a gradient in the number of e-folds of inflation. Thus, the scale of the cutoff in the 
spectrum varies spatially according to Eq. (fl~6|) . but we will approximate the spectrum of the 
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Figure 8: Averaged theta function spectra for the northern hemisphere for differing amounts of 
gradients such that ANi ow = 0, 1, 2, 3, 4, oo. 



uniform fast-roll model of Fig. [2]by a simple theta function cutoff spectrum P$(k) = Q(k — k c ). 
For simplicity, we do not include a spectral index and we compare the resulting spectra to a 
scale invariant Harrison-Zel'dovich spectrum. 

In the hemisphere analysis in Sec. S]we first calculated CMB spectra and then performed 
an average over the two hemispheres. Here however, we will first average the primordial 
spectrum over the hemispheres and then calculate the CMB power spectra. The primordial 
spectra for the northern and southern hemispheres are 



dz P$(k, z) 



dz P$(k, z) 



Zhigh 



Clow / Q 

d( exp -ANi ow - — 

V yom 

( c 

d( exp I -ANi ow - — 

(high V {low 



P*(fc,C) 



(27) 



from which we can calculate the averaged CMB spectra Cf and Cf. As seen in Eq. (|16p 
today's physical scale of the cutoff in the spectrum varies exponentially such that the spectrum 
becomes 

c 



k — ko exp +AN t 



Clo 



(28) 



where the position of the cutoff in the northern hemisphere varies from fco to k\ ow = koe ANlow 
depending on the position (. Plugging this into Eq. ([2?]) we obtain an average spectrum 



P£(k) 



1 - k /k 



6(fc - ko) @{h ow -k) + 0(k - h ow ) 



1 - k /ki ow 

for the northern hemisphere which vanishes for k < ko and equals unity for k > k\ ow . 
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Figure 9: Averaged theta function spectra for the northern (dashed) and southern (solid) hemisphere 
for AN low = 4. 



In Fig. [8] we display the averaged spectrum of the northern hemisphere for different 
amounts of gradient, i.e. for different ANi ow , where the limit ANi ow = corresponds to no 
gradient at all and we recover a pure theta function spectrum, whereas the limit AN[ ow = oo 
corresponds an infinite amount of gradient within the hemisphere. As one would expect, 
the average spectrum of the hemisphere is dominated by the part of the gradient with most 
e-folds of inflation which will dominate and make the average power spectrum rather steep 
around k$. At k/ko = 10 the spectrum already reaches 0.9 in the limit of an infinite gradient, 
and since i ~ k in the CMB spectrum, we expect that the cutoff in the resulting averaged Ci 
spectrum is also rather steep with a range of significant suppression up to roughly I ~ 10. 

The averaged spectrum for the southern hemisphere can be calculated analogously, and 
for ANi^ = 4, we show both averaged spectra in Fig. [9j The averaged spectrum of the 
southern hemisphere is steeper than the spectrum of the northern hemisphere since it contains 
a much smaller piece of the initial gradient with a difference of only about 0.69 e-folds inside 
the southern hemisphere as compared to ANi ow = 4 e-folds of difference in the northern 
hemisphere. 

As the next step we numerically calculate the CMB spectrum 0$, where for this toy model 
we will only account for the Sachs Wolfe effect [29] [3D] 



_ 4vr f°° dk . 2 f2k\ rtn . , Qn > 

ignoring all other effects such as the integrated Sachs- Wolfe effect. These Sachs- Wolfe spectra 
for the northern and the southern hemispheres are divided by the Sachs- Wolfe spectrum of a 
scale invariant Harrison Zel'dovich spectrum to find the suppression functions F^ and Ff for 
both hemispheres. For that we have to specify the scale Icq. We choose ko = clqHq such that 
the CMB spectrum of the northern hemisphere (z = . . . Z[ ow ) will exhibit a suppression of 
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Figure 10: Suppression functions resulting from an averaged theta function spectra for the northern 
(dashed) and southern (solid) hemisphere for ANi ow = 4. Note that the underlying calculation 
neglected the integrated Sachs- Wolfe effect. 

power at largest scales, whereas the CMB spectrum of the southern hemisphere z = Zhigh • • • 
will exhibit (almost) no large scale suppression. The resulting suppression functions and 

are shown in Fig. [TU1 and as one would expect, their shapes do not differ much from 
the primordial averaged power spectra in Fig. [9j At small scales or high £, both spectra 
match onto the CMB spectrum from a scale invariant primordial spectrum, and for observable 
multipoles £ > 2 the southern spectrum doesn't significantly vary from a featureless spectrum. 
The spectrum of the northern hemisphere however is significantly suppressed at low £, but 
as expected, for £ > 12 the suppression becomes smaller than 10%. 

We can test our theta function cutoff toy model by calculating the x 2 of the measured 
CMB data points in Fig. [5] with respect to the averaged spectra of the two hemispheres. 
The resulting x 2 f° r the northern hemisphere decreases by 2.6 in comparison to the isotropic 
ACDM model whereas x 2 doesn't change significantly for the southern hemisphere. Unlike 
in Sec. d] where we fitted the parameters and obtained Ax 2 = —3.1 for the best fit, here we 
simply chose plausible parameters ko and AN[ ow . 

For this toy model, we have only considered the Sachs- Wolfe contributions to the CMB 
anisotropies, ignoring the integrated Sachs- Wolfe effect etc. Nevertheless, when we compare 
the hemisphere suppression functions in Figs. [7|and[T0l we note that neither the precise form 
of the cutoff nor the inclusion of the integrated Sachs- Wolfe effect seem to change the main 
features. At the very lowest multipoles £ < 5 one can see in Fig. [7] the effect of the integrated 
Sachs- Wolfe effect resulting in a plateau instead of a cutoff extending down to zero. The main 
features however are determined by the geometric effect due to a gradient in the number of 
e-folds and by the steepness the cutoff feature in the spectrum. 

Before proceeding to study the next toy model, we will briefly comment on the integrated 
Sachs- Wolfe effect which we omitted in the analysis of the above toy model. As mentioned 
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Figure 11: Suppression functions from fully numerical CMB spectra resulting from uniform primor- 
dial theta function spectra with the cutoff position varied linearly in k. The plateau to the left of 
each cutoff arises due to the integrated Sachs- Wolfe effect, and its height depends on the location of 
the cutoff. 

earlier, the integrated Sachs- Wolfe effect gives contributions to Ci from a wider range of scales 
than the Sachs- Wolfe effect does. Therefore at scales larger than the scale associated with 
the cutoff, the Sachs- Wolfe effect contributes almost nothing to the CMB power spectrum 
but the integrated Sachs- Wolfe effect can still yield some power in that region. That gives the 
plateaus or offsets in the CMB spectrum. To illustrate these plateaus, it is most transparent 
to show the CMB spectrum resulting from a uniform theta function primordial spectrum, as 
seen in Fig. [11] for various cutoff scales. We observe that the height of the offset depends on 
the scale of the cutoff where smaller cutoff scales result in a lower offset. 

At the other extreme from a step function cutoff is one for which the cutoff can be made 
gradual. This can be modeled by use of the exponential cutoff spectrum introduced in Eq. 
(|10p where we allow for the exponent a to take on different values. The parameter a governs 
the steepness of the cutoff in the spectrum. For low values of a we obtain spectra which 
exhibit a slowly varying cutoff rather than the steep cutoff for a = 3.35. This latter value 
mimics the the CPKL fast-roll spectrum, and larger values than this would approach the step 
function. Here we explore the softer cutoff provided by small values of a. As we did for the 
theta function model, we keep the aspect of a gradient in the number of e-folds of inflation. 
The resulting scenario may be relevant for other models, with the basic features of an IR 
cutoff of the perturbation spectrum with the scale of the cutoff varying spatially as dictated 
by a gradient in the number of e-folds of inflation. 

As we did in the case of the theta function spectrum analysis, we do not include a spectral 
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Figure 12: Suppression functions resulting from averaged exponential spectra for the northern 
(dashed) and southern (solid) hemisphere for a = 0.5, £o = 2 and ANi ow — 2. Note that the 
underlying calculation neglected the integrated Sachs- Wolfe effect. 



index so that Eq. ()10p simplifies to 

P*(Jfe) = 1 - exp [ - (A;/A; c ) c 



(31) 



We find that when one numerically calculates the Sachs- Wolfe spectrum of such an expo- 
nential cutoff spectrum via Eq. (j30j) . the resulting Sachs- Wolfe spectrum is approximated 
extremely well by an exponential cutoff spectrum for £{£ + 1) C( (where the exponential a for 
the Sachs-Wolfe spectrum is a bit smaller than a in P$(k)). Ignoring again the integrated 
Sachs- Wolfe effect, we can therefore use 



£(£ + 1) C e = 1 - exp [ - (£/£ c f] 
where now the position of the cutoff in £ space varies spatially as 
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To find the averaged spectra for the two hemispheres, we integrate over the position dependent 
spectrum analogously to Eq. ([27|) 

d( exp ( -ANi ow 



Clc 



c 



Chigh 



d( exp ( -AN low 



C 



do 



Q(C) 



(35) 



20 



For the parameters a = 0.5, Iq = 2 and ANi ow = 2, the suppression functions FP for 
D = N,S are displayed in Fig. [I2j From this we see that it is possible to have a suppression 
factor that extends out to larger values of i, and which creates an asymmetry on the sky. 
Treated simply as a suppression of the standard WMAP spectrum, this modification will 
not lead to an improved fit because the suppression of the power spectrum occurs in both 
hemispheres. It is possible that a better fit could be obtained if one adjusts the magnitude of 
the power spectrum, but this would involve adjustments of all the astrophysical parameters 
in order to not upset the agreement obtained at smaller angular scales. 
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